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Abstract 

On the viewpoint of corpuscular model an electromagnetic radiation can be regarded as a system 
composed of photons with different energies and momenta, which provides us a method being 
different from the Maxwell wave theory to describe the interaction of electromagnetic waves with 
plasmas. In this paper the evolution behavior of a single photon and the collective effect of a 
photon system in plasma waves are uniformly described in the frame of photon dynamics. In a 
small-amplitude plasma wave the modulation of photon dynamical behavior by the plasma wave 
can be treated as perturbation, and the photon acceleration effect and photon Landau damping are 
investigated in the linear theory. In a plasma wave with arbitrary amplitude the photon evolution 
trajectories in phase space and coordinate space are obtained by solving the dynamical equations, 
and the trapping condition and possibility of photons in the given plasma wave are also discussed. 
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I. INTRODUCTION 



The interaction of intense fields with plasmas is an important hotspot issue in plasma 
physics, which is involved in many area, for example the excitation of wakefield, laser-plasma 
acceleration, fast ignition fusion[l-7], etc. A laser pulse can excite a plasma wave via the 
effect of ponderomotive force when propagating through a plasma, and motions of electrons 
and photons in plasmas are also modulated by this plasma wave. In the interaction process 
the laser-plasma system can be described by the Maxwell-fluid theory which has achieved 
huge success in the last two decades [8- 16]. 

On the viewpoint of corpuscular theory an electromagnetic field interacting with plasmas 
can be treated as a nonequilibrium photon system governed by the kinetic equation, which 
is called as photon kinetic theory [17-19]. Based on the photon kinetic theory the dynamical 
behavior of a single photon in a background plasma can be extended to the evolution of 
the entire photon system. As a special photon system, laser pulse has a lot of advantages 
compared with other electromagnetic radiations, such as good directivity, high intensity and 
narrow spectrum width, etc. We can treat the laser pulse as an ordered photon current 
described by the kinetic theory. The connection between the photon system and a single 
photon is the photon number distribution in phase space which is essentially the statistical 
weight of photons with various modes in the entire photon system. The evolution of an 
electromagnetic radiation in plasmas can be derived by making weighted average to all 
photon behaviors. Due to the photon kinetic theory being a very effective method for 
describing the propagation properties of electromagnetic fields in plasmas, the research on 
single photon dynamics is worth trying. 

There are some similarities between photon and electron behaviors in plasma waves. For 
instance, if the photon velocity matches with the phase velocity of the plasma wave, the 
photon can be accelerated effectively by the plasma wave[20-23] and the photon system 
can produce Landau coUisionless damping to the plasma wave[24-26], etc., and so do the 
electrons [27-31]. However the differences are obvious. The photon is not influenced by the 
electrostatic fleld of the plasma wave because of the electric neutrality, whereas the electron 
is. The propagation behavior of photons in plasmas is related to the refractive index, i.e., 
electron number density of the plasma. The refractive index of a plasma can be equivalent 
to a special potential fleld acting on photons, and photons will get some equivalent forces 
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when moving in this potential field just as electrons are affected by the electrostatic forces 
in plasma waves. In this sense, the dynamical evolution of a single photon in a plasma is 
determined by the electron number density of the plasma. In inhomogeneous plasma the 
equivalent forces acting on photons depend on their positions in the plasma, and the photons 
with different modes are also influenced by different forces owing to dispersion effects. 

The purpose of this paper is to investigate the propagation and evolution behaviors of 
photons in a background plasma wave by means of the photon dynamical theory. The photon 
dynamical equations are expressed in Hamiltonian formulation based on the approximation 
of geometrical optics [22], and the photon frequency is regarded as Hamiltonian for a photon. 
Based on the dynamical theory we describe some phenomena in the interaction process of 
photons with plasma wave in a unified frame. If the amplitude of a plasma wave is small 
enough, the modulation of the photon dynamical behavior by the plasma wave can be treated 
as perturbation. In Sec. 2 the photon dynamical equations are treated by the perturbation 
theory, and the photon acceleration effect, along with the collective coUisionless damping 
originated from the laser pulse, is analyzed to the second order. In Sec. 3 the dynamical 
evolution of a single photon in a plasma wave with arbitrary amplitude is investigated. The 
photon evolution trajectory in phase space is derived by solving photon dynamical equations. 
The trapping condition and possibility of photons in a given plasma wave are also analysed, 
and the motion equations in coordinate space for both trapped and untrapped photons are 
analytically obtained. At last we give a discussion in Sec. 4. 

II. PHOTON BEHAVIOR IN SMALL- AMPLITUDE PLASMA WAVE 

Let us consider one dimensional model for simplicity, and suppose that the electromag- 
netic wave propagates along x axis. In this paper we assume that the intensity of the 
electromagnetic wave is not strong enough to change the configuration of the plasma wave 
prominently so that the reaction from the electromagnetic wave (or the photon system) to 
the plasma wave can be ignored. When an electromagnetic wave propagates in a plasma 
wave, the disturbed electron number density modulates the frequency ui and wave vector k 
of the electromagnetic wave to produce frequency shifting, which can be described by the 
Maxwell wave theory. Nevertheless, we will discuss this problem in an alternative way. By 
considering the electromagnetic wave as a system of photon gas, let us discuss the kinetic 
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evolution of a single photon in a plasma wave. If the amplitude of a plasma wave is small, 
5n <^ no, where 5n is the disturbed electron number density and no is the ambient plasma 
density, the photon kinetics can be expressed by perturbation expansion: 

X ^ Xq + Vot + Xi + X2 + , (1) 

k = ko + h + k2 + , (2) 

U ^ Uo + COl + UJ2 + , (3) 

V = V0 + V1 + V2 + , (4) 

where x and v are the position and velocity of the photon, ou and k are regarded as the 
energy and momentum of the photon. The subscript denotes the initial quantity, and 1 
and 2 denote the first- and second-order corrections respectively, and Vq — c^ko/too- 

Assuming that the configuration of the plasma wave can be expressed by a function of 
/(x, t) = f{kpix — Upit), where kpi and Upi are the wave number and frequency of the plasma 
wave. The Hamiltonian of the photon is proportional to its frequency cu as determined by 
the dispersion relation of the electromagnetic wave in the plasma 

H = huj = h^c^k^ + uj^, (5) 

where oUp is the local plasma frequency determined by the disturbed electron number density 
in the plasma wave as oUp = oUpoi^ + f{x,t)). The photon evolution in a plasma wave is 
governed by the dynamical equations in Hamiltonian formulation 

dx OH 

dk _ _dH_ 

For simplicity, we suppose the plasma wave has the shape of cosine, / = K,cos{kpix — 
ujpit), where k, = 6n/no is the amplitude of the plasma wave and a small quantity k ^ 1. 
Dynamical equation (7) gives 

dk ul^kpiK . 

— = sin[kpix - ujpit). (8) 

Substituting Eqs.(l)-(3) into Eq.(8) the first- and second-order equations for the photon 
momentum are derived as 

dki ujpQkpiK . 

'dF ^ 2uJo ~ ' 
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dk2 
It 



2uJo 



kpiXi cos {kpiXo — ilt) sin {kpiXo — Hi) 

cos {kpiXo — fit) sin {kpiXo — fit) 



2uj, 







(10) 



in which fl — oUpi — kpiVQ. The integral of Eq.(9) gives a first-order correction of the photon 
momentum in the small-amphtude plasma wave 
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^pokpiK 

2ujQfl 



[cos {kpiXo — fit) — cos (kpiXo)] 



(11) 



In plasma waves the photon frequency is modulated by the electron number density, 
which is known as the photon acceleration effect. By using of the dispersion relation of an 
electromagnetic field in plasmas we can expand the photon frequency to second-order as 
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+ 
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cos^ {kpiXo — fit) — 



COS {kpiXo — fit) 



(12) 



4:C^kQUjQ ^ koUjQ 

where cuq and /cq are the initial photon frequency and momentum respectively. The frequency 
shift of photon in the plasma wave can be analyzed by using of Eq.(12). Defining the 
frequency shift at time t as Acu — uj{t) — uj{0), where uj{0) is the initial photon frequency, 
and using the first-order equation of photon momentum (11), it is given that 



'^pO'^plI^ 



[cos {kpiXQ + kpiAvt) - cos {kpiXo)] 



(13) 



2uJoAvkpi 

when calculated to the linear term, where the parameter Av = Vq — f^, and = Upi/kpi 
denotes the phase velocity of the plasma wave. If the photon velocity vq is very close to 
the phase velocity of the plasma wave, i.e., in the limit of Av 0, we get the maximum 
frequency shift 
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kpitujpQV^K, 
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sin (kpiXo) , 



(14) 



which is consistent with the earlier publications [21]. 

On the viewpoint of corpuscular property an electromagnetic radiation can be regarded as 
a particle system composed of photons with various modes. The evolution behavior of entire 
electromagnetic radiation in plasmas can be derived by making statistical summation of all 
photon behaviors in theory. Now we will discuss energy exchanging of the electromagnetic 



wave with the plasma wave based on the photon dynamics, which is a kind of collective 
effect of photon system known as photon Landau damping. 

Taking the derivative of Eq.(12) with respect to time t, the average rate of energy change 
of the photon to second-order in plasma wave is given by 

\dt / ujo \ uo J ujo 

{{xi sin [kpiXo — Qt)) — Q{xi cos {kpiXo — ^t))) 



2uo 



(^{ki cos (kpiXo — ^t) + Q{ki sin (kpiXo — ^t))^ . (15) 



2ul 

Here the dot stands for the derivative with respect to the time, and the bracket (• • •) denotes 
the average over initial position xq in one wavelength range of the plasma wave. The first- 
order correction of photon velocity can be obtained from dynamical equation Eq.(6) as 

C^ul^K (ul^kpi \ C^^pokpiK 

^ ^ 2ul \ u Q. ~ J ^ ~ ^ 2LJ^n ^^^^^pi^o)- (16) 

The integral of Eq.(16) gives 

^ 2ijjHl \^ ~ ^Vt ] ^^^^ ~ ~ (^pi^o)] 2a;° Q ^ ihi^o) ■ (1^) 

Substituting Eqs.(9)-(11), (16) and (17) into Eq.(15) we get 

^ y^^f^-^' s ( 1 _ (18) 

\dt I dn\n ^ y n ^ ^ ^ ^ 

In order to extend the rate of energy change of the photon to the whole electromagnetic 
field we need to computer the sum of all photons in the field. Supposing that the initial pho- 
ton number distribution is N{u!q), the average rate of energy change of the electromagnetic 
field is given by 

I (^)^ (^o) dcoo. (19) 

Eq.(18) indicates that (u) is oscillatory and damped rapidly, excepting at the resonance 
singularity = 0. Thus the main contribution in Eq.(19) comes from the area nearby the 
singularity, which can be worked out by considering the long time limit and expressed in 
terms of a delta function 

lim ^ sin (nt) = ttS (Q) = ^^^^5 (uoo - uj^) , (20) 



where a;<^ is resonance frequency of the photon with the plasma wave as defined by vo{ujq — 
— v^. Since d/dO. — — (wq^o) / i^po^pi) d/duo, we obtain 



t^'^\dt/ 8 duo ■ ^^^^ 

Inserting Eq.(21) into Eq.(19) the average rate of energy change of the electromagnetic field 
in plasma wave is expressed by 



TTOxIpok^CUpiK^ dN (cjo) 



(22) 
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Eq.(22) indicates a mechanism of energy exchange of the electromagnetic wave with plasma 
wave. The frequency a;^ can be regarded as the resonance frequency because the photon with 
this frequency has the velocity matching with the phase velocity of the plasma wave. Prom 
Eq.(22) we find that the energy exchange is determined by the derivative of the distribution 
function of photon number near a;^, and if the slope of this distribution function is steeper 
the energy exchange is more remarkable. That just means only the photons having the 
frequencies close to the resonance frequency are involved in this process of energy exchange 
with the plasma wave. Prom the preceding discussion it is clear that the energy exchange is 
not related to coUisions between photons and electrons. Thus the mechanism of the energy 
exchange here is a collisionless resonance effect between photon system and the plasma wave, 
namely the photon Landau damping[24,25]. Photon Landau damping has some properties in 
common with electron Landau damping: Both are collective velocity resonance phenomena 
with plasma waves. However the differences between them are obvious. The energy exchange 
in the process of electron Landau damping is achieved through doing work to electrons by the 
electrostatic field. But photons are electrically neutral and not influenced by the electrostatic 
force. In fact the energy exchange in photon Landau damping is achieved through the effect 
of photon frequency shift in plasma wave that is related to the fluctuation of the electron 
number density. 

The plasma wave is described by the Poisson equation 

= — ^ = _ 23 
Eq ax 

where and E are the scalar potential and strength of the longitudinal electric fleld in 
plasma wave. The solution of Eq.(23) is 

E — — - — sin {kpix — uipit) — E sin {kpix — uipit) . (24) 
s^okpi 



Then the amphtude of plasma wave k can be written as k — {eokpiE)/{eno), and Eq.(22) is 
given by 

\'^) ^zzz — ■ l^^J 



0)0 =W0 



Smno duo 

If we only consider the effect of photon Landau damping between photon system and the 
plasma wave and neglect other interaction, the damping rate 7 can be worked out by means 

of the energy conservation 



-SoPj + h{cu) = 0. 

Substituting Eq.(25) into (26) the damping rate is obtained as 

nchk^ojpikpi dN {uq) 
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(26) 



(27) 



Amrio dujo 

We consider the electromagnetic wave as a time-Gaussian laser pulse, the initial photon 
number distribution in frequency domain space is expressed by 



exp 



2 „2 



(wo - u;c) s 



(28) 



where Nq is a normalization constant depending on the intensity of the initial laser pulse 
and ujc is the center frequency and s is the duration time. Substituting Eq.(28) into (27) 
the damping rate of time-Gaussian laser pulse to the plasma wave is derived as 
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TT hrUpikpiUJcNo ck. 
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exp 



2 „2 



iO^) S 



(29) 



In Ref.26 the same damping rate is obtained by using the kinetic theory. If the phase velocity 
of the plasma wave satisfies the condition < Vc— (cujp) / (^c*^) , where Vc — c^l — Up/u'^, 
we have the damping rate 7 > for the time-Gaussian laser pulse. 



III. PHOTON EVOLUTION IN PLASMA WAVE WITH ARBITRARY AMPLI- 
TUDE 

Now we will discuss the situation of the plasma wave with arbitrary amplitude. It is 
more convenient to perform a coordinate transformation from the laboratory frame to a 
new frame: ^ — x — V(ft, t — t. In this new frame the shape of plasma wave has a simpler 
expression as / = f{kpi^), and the photon dynamical equations are derived from Eq.(6) and 
(7) as 

^^v-v^^ c^l-^ - v^, (30) 
8 



duj _ ujloV4,df{kpiC) , . 

dr~ 2uj ■ ^ ' 

In Eq.(31) we have canceled the photon momentum k by using of the dispersion relation of 
electromagnetic wave in plasmas. Combining Eqs.(30) and (31) the time r can be canceled 
and the evolution equation of photon in ^ — a; space is given by 

df ^1 - (1 + /)/VF - 

where ^^ — v^fy/c and W — u'^ /uj^q. Here we also call the ^ — u space as the phase space for 
convenience. The solution of Eq.(32) is 

y/^{yW-P^yJW-f-l)^l, (33) 

where a is a constant determined uniquely by the initial condition 

a = —— ^ (34) 

(\/Wo-/3^VW^o-/o-l) 

In Eq.(34) Wq depends on the initial photon frequency and /o is determined by the initial 
position of the photon injected into the plasma wave as /o — f{kpi^o). Eq.(33) indicates that 
(^\/W — (3^\/W — f — 1^ is an invariant as long as the initial conditions are given. From 
Eq.(33) the photon frequency is obtained as 




l±Pjl--o{l + f)\, (35) 



UJpO y/a 

which describes the cvohition trajectory of a photon in phase space. 

In Eq.(35) the notations "±" denotes the two different branches of the trajectory in 
phase space, which represents the evolutions of photons with two different kinds of initial 
conditions. The sign "+" stands for the upper branch and "— " denotes the lower branch. 
We also suppose the plasma wave has the shape of cosine, /(O = i^cos^kpi^), for simplicity. 
The photon evolution trajectory in phase space expressed by Eq.(35) is ilhistratcd in Fig.l, 
and the similar result was numerically obtained in Ref.22. According to the different initial 
conditions the evolution trajectory can be divided into two types: The first one is displayed 
as the dotted line in Fig.l, where the upper and lower branches intersect with each other 
to form a closed trajectory. The photons evolving in this trajectory will be confined to one 
wavelength area of the plasma wave and oscillate between the upper and lower branches. 
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FIG. 1: The photon evolution trajectory in phase space. The dotted hne denotes the evolution 
trajectory of trapped photon which is confined to one wavelength area of the plasma wave and 
oscillates between upper and lower branches. The dashed line denotes the evolution trajectory 
of untrapped photon, and these photons can propagate through different periods of plasma wave. 
The solid line is the separatrix between the trapped and untrapped photon trajectories. 

This type of photon is known as the trapped photon by the plasma wave. The second 
type is displayed as the dashed hne, where the upper and lower branches keep separated 
consistently. In this situation the photon evolves only on one branch and propagates through 
different periods of the plasma wave, so it is obviously not trapped by the plasma wave. If 
the photon is on the upper branch of the trajectory at the initial time, it will always evolve 
on the upper branch with the motion velocity faster than the phase velocity of the plasma 
wave. Whereas the photon will evolve on the lower branch and propagate backward with 
respect to the plasma wave if its initial position is on the lower branch in phase space. In 
Fig.l the solid line denotes the separatrix between trajectories of trapped and untrapped 
photons. Whether the photon can be trapped is determined by initial conditions. We have 
known that the trapped photon has the closed evolution trajectory, that is to say the factor 
under the square root sign in Eq.(35) 

rv 

1- ^{1 + Kcos{kpiO) =0 (36) 

70 

is solvable. Then we get the trapped condition for photons in plasma wave as 

<a< (37) 
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All the photons with initial conditions satisfying inequality (37) can be trapped by the 
plasma wave. The trapped photons have the velocity matching with the phase velocity of 
plasma wave, and can exchange energy effectively with the plasma wave. 

Let us discuss the motion equation of photons in coordinate space. The integral of 
dynamical equation (30) gives 

/ I = CT. (38 

^i-(i + /(e'))/W-(e') - P<t> 

Supposing that the plasma wave has the shape of cosine just as we discussed before and 

inserting Eqs.(33) and (35) into Eq.(38), the motion equation is given by 

r'^pi^ dx 



The sign "+" in Eq.(39) denotes the motion trajectory of the photon evolving on the upper 
branch, and "— '' represents the trajectory of the photon evolving on the lower branch. 

We will analyze the motion equation of trapped photon first. The maximum and minimum 
frequencies of the trapped photon in plasma wave can be obtained from Eq.(35) and given 

by 

u^a. = ^ {l + /3<^^1 - ;|(1 - «;) j , u^^ = ^ {l - /3^^1 - -|(1 - «:) j . (40) 

The initial conditions of the trapped photon satisfy inequality (37), which gives 1 — a/j"^ — 
ckk/j'^ < 0, and then the motion equation (39) is simplified into 

with 

" ( WtI) + W) - 1 ^ ^ 

The integral function in Eq.(41) can be expressed by the eUiptic integral of the first kind 
F{kpi^/2, k). For simplicity we only discuss the motion equation in the first period, the 
evolution after the first period can be worked out analogously. Assuming that the photon is 
on the upper branch of the trajectory in phase space at the initial time, the motion trajectory 
before the photon going into the lower branch is derived by using Eq.(41) as 

CT = P^lliC - eo) - ^^J^===lm [F{k,ii/2, k) - F{k,iio/2, k)] . (43) 
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And then, the photon reaches the lower branch and the motion equation is given by 
CT = /3^7?(e - Co) - , i^'l » _ ^ Im [2F{k,iW^, k) - F{kM2, k) - F{k,iC/2, k)] 



where determined by 



1 - 



;i + KCOs{kpi^rn)) = 



(44) 



(45) 



is the maximum position the photon can reach. At last the photon returns to the upper 
branch and finishes the evolution in a completed cycle. If the photon is on the lower branch 

of the trajectory in phase space at the initial time, we can analyze the motion equation by 
using the similar method, and the motion trajectory of this kind of photon is illustrated in 
Fig.2. It is shown that the trajectory is oscillatory, and the photon is confined to the area 





_J , , , , 1 , , , , 1 , , , , 1 , , , , cTk^ 

500 1000 1500 

FIG. 2: The motion trajectory of trapped photon in coordinate space. At the initial time the 
photon is on the lower branch of the trajectory in phase space. This kind of photon is confined to 
the area shorter than one wavelength of the plasma wave and oscillates between upper and lower 
branches, so the motion trajectory in coordinate space is also oscillatory as displayed in the figure. 

shorter than one wavelength of the plasma wave. The velocity of the photon is close to the 
phase velocity of the plasma wave. 

The situation of the untrapped photon is somewhat simpler. If the photon is on the 
upper branch of the trajectory in phase space at the beginning, it always propagate on the 
upper branch, and the maximum and minimum frequencies of the photon in plasma wave 
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are given by 



+ + (46) 



a 



Inequality (37) is not satisfied now and we have 1 — q;/7^ — olk/^"^ > 0, so the motion 
equation (39) gives 



27, 



cr = - 6) + T-^=^—Re [F(k,i^/2, ik) - F(k,i^o/2, ik)] 



with 



A;2 = 



(2««:/7') 



> 0. 



(47) 



(48) 



1 - («/7|) - («/^/7|) 

Eq.(47) is the motion equation of the untrapped photon with the initial position on the 
upper branch. In Fig. 3 the motion trajectory of this type of photon in coordinate space 
is displayed. The photon velocity is always faster than the phase velocity of the plasma 
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FIG. 3: The motion trajectory of trapped photon with the initial position on the upper branch in 
phase space. In this situation the photon always propagate with the velocity faster than the phase 
velocity of the plasma wave, so it travels forward with respect to the plasma wave. 

wave, so the photon propagates forward with respect to the plasma wave and can enter into 
different periods of the plasma wave. If the photon is on the lower branch of the trajectory 
in phase space at the initial time, it will always propagate on the lower branch, and the 
maximum and minimum frequencies of the photon in plasma wave are obtained by using 
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Eq.(35) as 



7| 



(49) 



The motion equation of this kind of photon is given by 

27^ 



Re [F(A;p,e/2, ik) - F{kM2, ik)] . (50) 



We illustrate the motion trajectory in Fig.4. The velocity of this type of photon is always 
slower than the phase velocity of the plasma wave, and the photon moves backward with 
respect to the plasma wave. In Fig.5 we show the motion trajectory of the photon evolving 




FIG. 4: The motion trajectory of trapped photon with the initial position on the lower branch in 
phase space. This kind of photon always moves with the velocity slower than the phase velocity of 
the plasma wave, so it travels backward with respect to the plasma wave. 



on the separatrix in one complete cycle. The photon is just restricted to one wavelength 
area of the plasma wave and oscillates between two neighbor peaks of the plasma wave. 

At last let us discuss a simple case in the limit of v^p — > c. The evolution equation of a 
photon in phase space, Eq.(33), is simplified as 



(51) 



in this situation, which gives 



{aKCOs{kpi^) + a + 1) 



(52) 
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FIG. 5: The motion trajectory of the photon evolving on the separatrix in one complete cycle. 

Considering the region of underdense plasma u ^ oOpQ, i.e., W ^ 1, both 1/a and {l + f)/W 
are small quantities, and the motion equation (38) can be expanded in series and calculated 
to the first order, which gives 

aKsm{kpi^) + (a + 2)kpi^ + Icrkpi = aKsm{kpi^o) + (a + 2)kpi^Q. (53) 

Apparently, in this limit the photon cannot be trapped by the plasma wave because its 
velocity is slower than c all the time, and it will propagate backward with respect to the 
plasma wave. 

IV. CONCLUSION AND DISCUSSION 

In this paper some dynamical behaviors of photons and collective effect of photon system 
in plasma wave are described in the unified frame of the photon dynamical theory. We treat 
the plasma wave as a special background for photons, and the refractive index of the plasma 
is equivalent to an external potential field which determines the dynamical evolution of 
photons. The dynamical equations of photons arc constructed in Hamiltonian formulation in 
the approximation of geometrical optics. We solve the dynamical equations and analyze some 
motion behaviors of photons in plasma wave, including photon acceleration and trapping 
effects, the evolution trajectory in phase space and photon Landau damping to the plasma 
wave. In fact the photon trapping effect in plasma wave is related to the photon acceleration. 
Only when the photon velocity is close to the phase velocity of plasma wave, the trapping 
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phenomenon occurs, which is analogous to the electron trapping effect in a plasma wave. 
Whether a photon can be accelerated to a velocity that matches with the phase velocity of 
a given plasma wave is determined by the initial conditions. The trapping condition and 
possibility of photons in a given plasma wave are analysed in detail in this paper. In a 
small-amplitude plasma wave, the evolution behavior of a single photon is extended to the 
entire electromagnetic field to calculate the coUisionless Landau damping effect originated 
from the electromagnetic field. When the plasma wave gets a relativistic phase velocity, for 
example a plasma wave driven by a laser pulse, some photons of the electromagnetic field 
can propagate with velocities very close to the phase velocity of the plasma wave and be 
trapped by the plasma wave, and consequently the velocity resonance with plasma wave 
occurs. The energy can be exchanged effectively between the electromagnetic field and the 
plasma wave in this resonance process. If we consider the electromagnetic field as a laser 
pulse, due to the narrow spectrum width most photons in the pulse may resonate with the 
plasma wave to produce the photon Landau damping effect. 

Theoretically, if we know the dynamical evolutions of photons in plasma waves, the 
evolution properties of the laser pulse can be obtained by making statistical summation of 
all photon behaviors, in which the keypoint is the weight factor of photons with different 
modes in entire pulse, i.e., photon number distribution. In a small-amplitude plasma wave 
the modulation on photon motion by the plasma wave is so weak that photon number 
distribution can be considered as a invariant, so we can use the initial pulse distribution as 
the statistical weight and keep it unchanged. However in a large- amplitude plasma wave 
the photon motion is significantly modulated and the photon number distribution evolves 
to obey the kinetic equation. Consequently, it is also an important problem to solve the 
photon kinetic equation in a large-amplitude plasma wave so as to obtain the time-dependent 
photon number distribution in phase space [32]. The photon dynamics plays an important 
role in the research of the interaction of electromagnetic field with plasma wave based on 
the corpuscular theory of light. By combining the kinetic theory of photon system and the 
single photon dynamics, the evolution of an electromagnetic field in a given plasma wave 
can be well- described in a different method from the Maxwell wave theory. 
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